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摘 Æ. 采用 2002—2016 年 GRACE (Gravity recovery and climate experiment) € 77 T Æ JPL-RLO6M 
数据 分 析 咸 海流 域 陆 地 水 储量 变化 (Terrestrial water storage change, TWSC) 时空 变化 特征 ,并 结合 
CRU TS4.03 气象 数据 \GLDAS-Noah 地 表 蒸 散发 数据 和 高 精度 土地 利用 数据 探究 气候 变化 与 人 类 活 
动 对 陆地 水 储量 的 影响 。 结 果 表 明 :(1) 2002 一 2016 年 成 海流 域 陆地 水 储量 变化 呈现 -3.20 mm-a” 
下 降 趋 势 , 春 、 夏 季 陆 地 水 储量 呈 有 蛋 余 态势 ,秋冬 季 呈 现 亏 损 状 态 ; 水 储量 变化 在 空间 上 表现 为 中 
部 和 东部 蛋 余 ,周边 亏损 的 特征 。(2) 2002 一 2016 年 威海 流域 降水 量 呈 -1.14 mm-a” 下降 趋势 ,地 


表 温 度 呈 0.11 Ca 上 升 趋势 ; 相 比 气温 ,水 储量 


i=! 


变化 与 降水 量 相 关 性 更 强 。(3) 2000—2015 年 ,成 


海流 域 耕 地 面积 小 幅 增 加 1.65x10* km ,水 域 面 积 减 少 ;农作物 耗 水 和 灌溉 需 水 的 增加 加 剧 了 威海 
流域 水 量 支出 ,威海 流域 燕 散 发 呈 21.63x10s ma :增加 态势 ,在 空间 上 与 陆地 水 储量 变化 的 相关 


系数 最 高 达 0.74, 是 影响 陆地 水 储量 变化 的 主 


EX —. 


X SE: 成 海流 域 ; 陆地 水 储量 ; 气候 变化 ; 土地 利用 变化 


文章 编号 : 


威海 曾 是 中 亚 最 大 的 内 陆 湖 ,水面 面积 曾 超过 


分 析 ,很 难 准确 评 佑 区域 水 资源 的 时 空 变 化 。 近 年 


6.70x10* km ,流域 面积 约 190x10' km。 由 于 远离 海 
TE ,其 水 资源 以 积 雪 和 冰川 融雪 补给 为 主 , 水 文系 
统 脆 弱 ”…, 对 气候 变化 的 响应 极为 敏感 汪 。20 世 纪 
60 年 代 以 来 , 阿 姆 河 与 锡 尔 河流 域 大 规模 的 水 土 资 
源 开 发 和 不 可 持续 的 灌溉 方式 ,导致 咸 海 湖面 快速 
EHR ,生态 环境 急剧 恶化 。 戌 海 危 机 ?已 成 为 20 世 
纪 人 类 活动 引起 的 最 严重 生态 环境 问题 之 一 ,受到 
国际 社会 的 广泛 关注 。 充 分 了 解 威海 流域 水 资源 
的 交换 过 程 和 赋 存 方式 , 探 明 区 域 水 资源 储量 及 其 
时 空 演 变 规律 ,对 成 海流 域 水 资源 管理 和 生态 保护 
有 重要 意义 。 

以 往 咸 海 流域 水 资源 变化 研究 多 基于 传统 气 
象 水 文 观测 资料 ,在 观测 站 点 稀少 且 观 测 不 连续 
的 条 件 下 ,对 湖泊 地 区 降水 量 蒸 发 .径流 等 进行 
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来 利用 多 源 遥 感 数 据 ,结合 地 基 观 测 或 建立 看 合 模 
型 水量 平衡 模型 ,可 提取 地 表 水 循环 要 素 的 时 空 
变化 信息 ,但 主要 是 对 水 域 面积 和 水 位 变化 ,以 
及 车 发 "一 等 特定 成 分 研究 ,而 对 陆地 水 循环 过 程 
中 河川 径流 、 地 下 水 等 要 系 缺 乏 适 用 性 。2002 年 发 
Se} 的 GRACE (Gravity recovery and climate experi- 
ment) 重 力 卫星 ,能够 可 以 连续 .快速 ,大 尺度 地 探测 
由 地 球 表 面 和 内 部 质量 迁移 造成 的 重力 变化 ,进而 
监测 地 表 水 储量 ”和 地 下 水 "变化 ,目前 已 经 成 
功 应 用 于 监测 全 球 “ 区域 ”和 流域 ”尺度 的 陆 
地 水 变化 。 基 于 GRACE 卫星 时 变 重 力 场 的 地 表 质 
量变 化 反 演 方法 按照 反 演 原理 的 不 同 主要 分 为 2 
类 : 球 谐 系数 反 演 方法 和 Mascon(Mass concentra- 
tion) 方 法 。 球 谐 系 数 法 原理 简单 .编程 实现 容易 ， 
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球 谐 系 数 法 反 演 结果 存在 严重 的 南北 方向 条 人 带 误 
差 ,必须 对 球 谐 系数 进行 滤波 处 理 ””。 基 于 点 质 
量 模型 原理 构建 卫星 间距 离 变 率 观测 值 和 地 表 区 
域 网 格 质 量变 化 的 直接 子 数 关系 式 ,建立 Mascon 方 
法 ;已 有 人 研究 表明 ,与 球 谐 系 数 法 反 演 结果 相 比 ， 
Mascon 模型 解 能 提高 反 演 结果 信 噪 比值 ,空间 分 辩 
率 更 高 ,保留 更 多 有 效 信号 ,与 GLDAS(Global land 
data assimilation system) 水 文 模型 结果 相关 性 更 好 二 ; 
^ [n] Mascon 模型 解 中 ,JPL(Jet propulsion laborato- 
ry) 发 布 的 Mascon 模 型 解 信 咯 比值 较 大 ,与 GLDAS 
水 文 模型 时 间 序列 相关 性 较 好 。 

因此 ,本 文选 取 中 亚 威海 流域 为 典型 研究 区 , 
采用 2002 年 1 月 至 2016 年 12 月 GRACE 陆地 水 储 
量 (Terrestrial water storage , TWS) Mascon 解数 据 , 结 
合 GLDAS-Noah 模型 的 地 表 蒸 散发 数据 以 及 CRU 
(Climate Research Unit) 发 布 的 降水 .气温 数据 ,分 析 
了 威海 流域 及 其 各 分 区 近 15 a 陆地 水 储量 变化 
(Terrestrial water storage change, TWSC) 趋势 ,并 探 
讨 了 气候 变化 与 人 类 活动 对 其 影响 。 研 究 结果 将 
有 助 于 从 宏观 上 了 解 中 亚 内 陆 湖 流域 水 资源 变化 
特征 ,为 水 资源 的 合理 分 配 与 管理 提供 依据 。 


1 研究 区 概况 


威海 流域 位 于 亚 欧 大 陆 腹 地 , 介 于 38"70' ~ 
62°20'E, 41°16'~46°52'N ZZ [B] ,总 面积 约 123x10; 
km (图 1)。 流 域 地 貌 格 局 清晰 ,变化 多 样 ,地 形 呈 
现 东 高 西 低 的 特点 。 西 部 和 西北 部 为 图 兰 平原 ,其 
间 分 布 卡 拉 库 姆 沙漠 和 克 孜 勒 库 姆 沙漠 ;东部 和 东 


35*E 60°E 65°E 70°E 75°E 


阿 姆 河 流域 上 游 m 锡 尔 河流 域 上 游 8 
加 阿 姆 河 流域 中 游 mm 锡 尔 河流 域 中 游 


图 例 
吕 威海 流域 em 阿 姆 河流 域 下 游 mm 锡 尔 河流 域 下 游 
55°E 60°E 65°E 70°E TPE 


图 1 威海 流域 
Fig. 1 The Aral Sea Basin 


南部 为 天 山 和 帕 米 尔 高 原 高 山区 ,是 主要 的 产 流 
De) 流域 具有 典型 的 大 陆 性 气候 特点 ,年 平均 
坡 , 超 过 2000 mm, 而 成 海 附 近 和 沙漠 区 年 平均 降雨 
量 不 足 100 mm**1。 

地 表 水 主要 由 阿 姆 河 和 锡 尔 河 两 大 水 系 构 
成 。 其 中 , 阿 姆 河 发 源 于 阿富汗 与 元 什 米尔 地 区 区 
界 处 兴 都 库 什 山脉 北 坡 维 略 夫 斯 基 冰 川 “ ,从 源头 
到 铁 尔 梅 效 灌区 划分 为 阿 姆 河流 域 上 游 ,以 东西 走 
回 的 高 大 山脉 为 主 , 中 下 游 分 区 以 Tujamuin 水 文 站 
点 为 界 , 下 游 为 地 势 平 坦 的 殉 漠 绿洲 ,流域 内 地 势 
落差 较 大 ,海拔 由 西向 东 逐 渐 升 高 ”。 锡 尔 河 发 源 
于 吉尔 吉 斯 境内 天 山 山 区 ,上 游 的 2 条 主要 文 流 纳 
伦 河 与 库 瓦 河 在 费 尔 干 纳 汇 合 后 称 为 锡 尔 河 ” ,从 
河源 到 费 尔 干 纳 盆地 出 口 为 锡 尔 河流 域 上 游 , 费 尔 
干 纳 盆 地 出 口 至 恰 尔 达 拉 水 库 为 中 游 , 恰 尔 达 拉 水 
库 至 河口 为 下 游 。 两 河流 域 中 下 游 地 区 适宜 的 气候 
条 件 和 丰富 的 土地 资源 ,非常 适合 发 展 灌溉 农业 “。 
20 世 纪 60 年 代 以 来 ,由 于 乌 效 别克 斯 坦 和 哈萨克 斯 
坦 两 国 大 规模 的 耕地 扩张 ,改变 了 威海 流域 的 水 量 
平衡 状态 ,水 资源 配置 的 变化 导致 两 河 注 入 咸 海水 
量 显著 减少 , 咸 海 湖泊 面积 急剧 萎缩 、 水 位 持续 降 
低 ,生态 环境 问题 愈演愈烈 。 


2 数据 与 方法 


2.1 数据 来 源 

GRACE 陆地 水 储量 数据 为 JPL 发 布 的 RLO6M 
数据 ,时 间 跨 度 为 2002 年 1 月 至 2016 年 12 月 共 180 
个 月 ,由 于 卫星 目 身 原因 和 测量 误差 等 因素 的 影 
啊 ,本 文 实际 获得 该 时 段 的 156 个 月 份 数 据 ,缺失 24 
个 月 份 数据 。 为 了 保证 数据 结果 的 准确 性 ,将 24 个 
缺失 月 份 的 数据 记 为 空 值 ,不 进行 插 补 ,空间 分 辨 
X H 0.5°x0.5° ,时 间 分 辩 率 为 1 个 月 (https://grace. 
jpl.nasa.gov/data/get-data/jpl_global_mascons/ ) 。 

GLDAS 陆 面 同化 数据 为 Noah 模 式 模拟 的 陆 表 
蒸 散发 数据 (https://disc.gsfc.nasa.govw) 产 品 ,空间 分 
辨 率 为 1"xl?" ,时 间 分 辩 率 均 为 3h, 下 载 数据 时 段 与 
GRACE 卫星 数据 一 致 。 

气象 数据 采用 英国 East Anglia 大 学 CRU 发 布 
的 CRU TS4.03 数据 (http://data.ceda.ac.uk/bade/cru/ 
data/cru_ts/) , F ZI] Ex 7y 2002 Æ 1 H Æ 2016 4 12 
月 ,空间 分 辨 率 为 0.5°x0.5° 的 气温 和 降水 月 数据 。 
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土地 利用 数据 采用 中 国 科 学 院 " 地 球 系统 科学 
数据 共享 平台 -新 疆 与 中 亚 科 学 数据 共享 平台 ”发 
布 的 2000、2015 年 中 亚 土 地 利用 数据 ,该 数据 集 参 
照 2000 年 TM 影 像 和 2015 年 OLI 影 像 , 采 用 人 工 目 
视 解 译 技术 进行 修改 和 确认 后 获得 。 解 译 结果 将 


式 中 : 设 定 有 x 和 ;两 组 要 素 的 n 对 观测 资料 xi、 
Mat xS Vi Yost us n EEA x 是 变量 x 的 平 
均值 ;y 是 变量 y 的 平均 值 。 当 7 为 正 值 时 ,表示 x 
与 y 之 间 正 相关 ,rr 为 负 值 时 ,表示 x 与 y 之 间 负 相 
Ko 一般 情况 下 ,lrl=1.0 表 示 完 全 相关 ;0.8<Irl<1.0 


人 研究 区 分 为 耕地 林地、 草地 城乡 居民 工矿 用 地 、 
水 域 . 未 利用 土地 6 个 一 级 类 。 

2.2 研究 方法 

2.2.1 GRACE 时 变 重 力 场 反 演 水 储量 变化 计算 方 
法 Watkins 57 FES AERA Re f 4551 个 质量 
均匀 分 布 的 3° 等 面积 球 帽 (MASCON) ,利用 JPL 重 力 
场 Level-1 数 据 , 通 过 加 权 最 小 二 乘 方法 得 到 每 个 球 
帽 的 等 效 水 柱 高 值 。 本 研究 基于 GRACE 卫星 时 变 重 
力 场 模型 数据 输出 为 经 纬 网 节点 数据 ,利用 ArcGIS 
平台 进行 区 域 统计 分 析 , 所 需要 的 工具 包括 : Make 
NetCDF Raster Layer, Zonal Statistic as Tabled, Band 
Collection Statistics , 4f HH Yat ERK fii e AE eA I, 38 
过 计算 得 到 各 年 月 份 流 域 水 储量 变化 统计 表 ™。 
2.2.2 时 间 序 列 趋势 分 析 方法 ”本 文采 用 线性 趋势 
法 描述 陆地 水 储量 、 降 水 量 、 气 温 、 燕 散 量 时 间 序 询 
的 趋势 性 变化 。 对 流域 内 CRU 降水 量 和 气温 数据 
以 及 GLDAS-Noha 1*x1° 蒸 散发 数据 进行 2002 一 
2016 年 际 变化 空间 趋势 分 析 。 长 期 变化 趋势 的 显 
著 性 检验 采用 Mann-Kendall 检 验 法 3。 

2.2.3 空间 相关 性 分 析 陆地 水 储量 变化 与 降水 量 
以 及 地 表 温 度 之 间 关 系 的 研究 采用 基于 栅 格 的 空 
间 相 关 性 分 析 方 法 ,通过 逐个 像 元 计算 陆地 水 储量 
与 降水 量 .地表 温度 等 之 间 的 相关 系数 来 表征 陆地 
水 储量 与 其 他 因子 的 相关 性 ,计算 公式 如 下 : 


n 


S [e — xy; - y) 
E : (1) 
Se. - x Y ty, =F) 


i=l i=l 


100r (a) 年 际 变化 
y=—0.2663x+14.107 


陆地 水 储量 变化 /mm 


-1005002 2004 2006 2008 2010 2012 2014 2016 
年 份 


陆地 水 储量 变化 /mm 


表示 高 度 相 关 ;0.5<lrl<0.8 表 示 显 著 相 关 ;0.3<lrl<0.5 
表示 低 相 关 ; 0.0<lrls0.3 表 示 微 弱 相 关 ;lrl=0.0 表 示 
完全 不 相关 。 


3 结果 与 分 析 


3.1 威海 流域 陆地 水 储量 时 空 变 

3.1.1 陆地 水 储量 的 年 际 和 年 内 变化 图 2a~b 展 
示 了 2002 年 1 月 至 2016 年 12 月 成 海流 域 陆地 水 储 
量 年 际 和 年 内 变化 。 从 年 际 变化 来 看 , 近 15 a 咸 海 
流域 水 储量 整体 以 -3.20 mm .a-' 速 率 呈 显著 下 降 趋 
势 (P<0.05) ; 最 大 盘 余 出 现在 2005 年 3 月 ,高 达 
81.05 mm, 最 大 亏损 出 现在 2008 年 10 月 ,为 -91.29 
mm( 图 2a)。 从 年 内 变化 来 看 ,陆地 水 储量 变化 的 
峰值 多 出 现在 4 月 或 5 月 。 春 季 (3 一 5 月 ) 和 秋季 
(9—11 H ) 陆 地 水 储量 变化 高 于 夏季 (6 一 8 月 ) 和 冬 
4 (12 H —IK4E 2 A) (Al 2b). 

3.1.2 陆地 水 储量 的 室 间 变化 ”2002 一 2016 年 威海 
流域 陆地 水 储量 呈现 中 部 腹地 和 东部 鱼 余 ,周边 亏 
损 的 特征 ;其 中 流域 北部 和 西部 水 量 亏 损 程度 大 于 
南部 地 区 ,特别 是 在 土库曼 斯 坦 境内 及 乌 效 别克 斯 
坦 境 内 的 威海 地 区 ,陆地 水 储量 的 亏损 更 为 显著 ， 
而 水 储量 的 和合 余 则 主要 出 现在 流域 东部 山区 以 及 
阿 姆 河 流域 中 游 的 灌区 (图 3)。 

为 了 解 咸 海流 域 陆 地 水 储量 的 时 空 变异 特征 ， 
将 咸 海 流域 划分 为 阿 姆 河 流域 和 锡 尔 河流 域 ,并 依 
据 流 域 地 貌 、 水 文 站 点 和 灌区 分 布 , 划 定 了 阿 姆 河 
流域 和 锡 尔 河流 域 的 上 中、 下 游 分 区 。 按 照 区 域 


(b) 年 内 变化 
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图 2 2002 一 2016 年 威海 流域 陆地 水 储量 年 际 和 年 内 变化 
Fig. 2 Interannual and intra-annual changes of TWSC in the Aral Sea Basin from 2002 to 2016 
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划分 提取 陆地 水 储量 的 变化 序列 , 结 末 表 明 : 阿 姆 
河流 域 和 锡 尔 河流 域 水 储量 均 呈 显著 下 降 趋势 (图 
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图 3 2002 一 2016 年 成 海流 域 陆地 水 储量 变化 空间 分 布 
Fig.3 Spatial distribution of TWSC in the Aral Sea Basin 
from 2002 to 2016 
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4a-b) ,变化 速率 分 别 为 -2.91 mma! 和 -5.18 mm-a” 
(CP<0.05,2008 年 后 P<0.01)。 除 阿 姆 河流 域 中 游 水 
储量 变化 呈 1.79 mm:a 小幅 增加 外 (图 4e) ,其 他 5 
个 分 区 陆地 水 储量 均 呈 减 小 趋势 ;其 中 阿 姆 河流 域 
下 游 减 小 幅度 最 大 ,变化 速率 达 -7.36 mm a (Al 
4g) , 阿 姆 河 流域 上 游 和 锡 尔 河流 域 上 游 分 别 以 
-7.28 mm.a -3.81 mm: a 速率 减 小 (图 4c~d)。 锡 
尔 河 流域 上 .中 .下游 陆地 水 储量 减 小 幅度 为 :上 游 
< 下 游 < 中 游 ( 图 4d、4f、4h)。 
3.2 环境 变化 对 威海 流域 陆地 水 储量 的 影响 
3.2.1 咸 海流 域 陆 地 水 储量 对 气候 要 素 的 敏感 性 
气候 变化 与 人 类 活动 是 影响 威海 流域 水 储量 变化 
的 主要 因素 。 本 研究 选取 降水 量 .气温 为 主要 气候 
子 , 分 析 其 与 陆地 水 储量 的 关系 。 
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图 4 威海 流域 及 其 分 区 陆地 水 储量 的 年 际 变化 
Fig.4 Interannual change of TWSC in the Aral Sea Basin and its sub-region 
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2002 一 2016 年 咸 海流 域 降水 量 整 体 呈 减 小 趋 
势 ( 图 5a) , 减 小 速率 约 为 -1.14 mm: a! ;减少 最 明显 
在 成 海流 域 东 部 (-8.22 mm * a^) ;降水 量 增加 最 明 
显 的 地 区 出 现在 锡 尔 河流 域 和 阿 姆 河 流域 的 下 游 ， 
增幅 可 达 2.67 mma (B| 5c) ; 2002—2016 4 Jk 
域 的 地 表 温 度量 上 升 趋势 ,上 升 速率 为 0.11 Cea! 
(图 $b) ,其 中 下 降 明 显 的 地 区 为 阿 姆 河 流域 下 游 中 
部 , 达 -0.03 Cea! ;地 表 温 度 最 大 增幅 在 阿 姆 河流 域 
南部 约 0.05 WC.a!'( 图 5d)。 咸 海流 域 陆 地 水 储量 和 
降水 量 的 相关 系数 整体 上 多 呈正 相关 关系 ,中 下 游 
灌区 和 南部 以 及 东部 山区 的 相关 系数 较 高 (图 5e); 
正 相 关 出 现在 流域 上 游 和 中 部 区 域 , 负 相 关 主 要 出 
现在 锡 尔 河流 域 下 游 与 阿 姆 河流 域 中 游 地 区 和 阿 
姆 河流 域 下 游 的 少 部 分 地 区 。 陆 地 水 储量 与 地 表 
温度 在 正 相 关 区 域 主要 集中 在 锡 尔 河流 域 中 下 游 
地 区 以 及 与 阿 姆 河 流域 中 下 游 交 汇 处 ,而 负 相 关 


(a) 年 降水 量 
350 y= —1.137x*253.24 


年 降水 量 /mm 


2002 2004 2006 2008 2010 2012 2014 2016 


年 份 
(c) 降水 量变 化 趋势 


则 主要 分 布 在 两 河流 域 的 上 游 地 区 以 及 阿 姆 河流 
域 的 南部 和 下 游 的 西南 部 (图 5f)。 总 的 来 说 ， 
2002 一 2016 年 威海 流域 降水 减少 导致 水 量 收入 减 
少 , 而 气温 升 高 可 引起 蒸 散发 增加 导致 水 量 文 出 增 
加 ,这 与 陆地 水 储量 的 整体 减 小 趋势 一 致 ;降水 量 
和 气温 在 一 定 程 度 上 影响 流域 的 陆地 水 储量 的 季 
节 性 变化 ,但 气温 与 陆地 水 储量 的 相关 系数 和 显著 
性 并 不 高 ; 相 比 气温 ,水 储量 与 降水 量 的 相关 性 
更 强 。 

3.2.2 土地 利用 变化 对 成 海流 域 陆 地 水 储量 的 影 
-& 咸 海 流域 5 个 国家 的 灌区 面积 从 1970 年 开始 
迅速 增长 ,1991 年 苏联 解体 后 ,各 国 水 资源 开发 利 
用 与 管理 体制 发 生 了 很 大 变化 ,耕地 开 屡 的 强度 和 
速度 较 前 期 而 言 都 大 幅 降 低 , 以 缓慢 增加 和 小 幅 减 
少 为 主 ,但 以 灌溉 为 主 的 农业 耗 水 量 依 然 很 大 “。 
图 6 为 2000 年 和 2015 年 威海 流域 土地 利用 窗 盖 变 
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图 5 威海 流域 气象 要 素 变 化 及 其 与 陆地 水 储量 变化 的 相关 关系 


Fig. 5 The variation of meteorological elements and the correlation with TWSC in Aral Sea Basin 
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图 6 2000 年 和 2015 年 威海 流域 土地 利用 图 盖 变化 
Fig. 6 Changes of land use cover in the Aral Sea Basin from 2000 to 2015 


化 图 。 与 2000 年 相 比 ,2015 年 威海 流域 耕地 、 林 地 、 
城乡 工矿 居民 用 地 面积 均 有 所 增加 ( 表 1) ,成 海流 
域 耕地 面积 增加 了 1.65x10* km?, 其 中 阿 姆 河流 域 耕 
地 增加 了 1.19x10' km , 锡 尔 河流 域 增 加 了 0.45x10” 
ko’; 草地 水域 面积 减少 ,分 别 减 少 了 1.52x10' km? 
All 1.56x10* km’. 

2002—2015 年 阿 姆 河流 域 耕 地 面积 增加 都 是 
主要 的 变化 ,特别 是 在 阿 姆 河 流域 中 游 地 区 ,2000 
年 至 2015 年 新 增 耕 地 面积 0.76x10' km , 占 阿 姆 河 
流域 耕地 增加 总 量 的 63.87%。 不 同 的 是 ,2015 年 锡 
尔 河 中 游 耕 地 面积 较 2000 年 小 幅 减少 47 km’, 5:55 
尔 河流 域 耕地 面积 变化 的 1.04% , 阿 姆 河 流域 各 分 
区 增加 面积 总 体 大 于 锡 尔 河流 域 (图 7a)。 此 外 ,水 


已 有 研究 表明 ,1960 一 2015 年 咸 海面 积 与 流域 
灌溉 面积 和 水 库容 量 旦 显著 的 负 相 关 关 系 中 。 土 地 
利用 覆盖 类 型 和 面积 的 变化 ,改变 了 水 资源 在 流域 
中 的 时 空 配置 ,从 而 改变 了 整个 流域 水 资源 的 平衡 
状态 。 最 直接 的 表现 就 是 蒸 散发 的 变化 ,基于 GL- 
DAS 蒸 散 数 据 的 分 析 表 明 , 2002—2016 年 威海 流域 
年 蒸 散 量 以 21.63x10: m? a^ 速率 增加 ; 阿 姆 河 流域 
和 锡 尔 河流 域 蒸 散 量 分 别 以 11.02x108 m?* a! ,8.36x 
105 na 的 速率 增加 。 伴 随 着 耕地 和 水 域 面 积 的 
增加 (图 7a) ,各 分 区 蒸 散 量 除 阿 姆 河流 域 下 游 以 
—0.26 10° m° -a 的 速率 下 降 之 外 ,其 他 5 个 分 区 蒸 
散 量 均 呈 不 同 幅度 增加 (图 7b~d)。 而 阿 姆 河流 域 
下 游 则 由 于 威海 萎缩 导致 水 面 蒸 发 持续 减 小 ,区 域 


域 面 积 的 变化 也 较为 显著 , 阿 姆 河流 域 由 于 南 成 海 
的 持续 萎缩 水 体面 积 减 小 ,而 锡 尔 河流 域 水 体面 积 
则 小 幅 增 加 ;其 中 , 锡 尔 河 下 游 分 区 水 域 面积 增加 
最 多 , 约 0.12x10’' km’, 


蒸 散 量 呈 下 降 趋 势 。 锡 尔 河流 域 上 下游 地 区 农 作 
物种 植 面积 增加 导致 蒸 散 量 增 大 ;而 中 洲 地 区 耕地 
和 水 域 面 积 变化 很 小 , 花 散 发 的 变化 可 能 更 多 地 受 
气候 变化 的 影响 ,升温 导致 燕 散 量 增加 。 


表 1 2000 年 和 201S 年 威海 流域 及 其 子 流 域 土地 利用 面积 变化 


Tab. 1 Changes in land use areas in the Aral Sea Basin and its sub-basin in 2000 and 2015 /10* km? 
n 威海 流域 阿 姆 河流 域 锡 尔 河流 域 
2000 年 2015 年 2000 年 2015 年 2000 年 2015 年 
耕地 12.91 14.56 7.94 9.13 4.98 5.43 
林地 2.89 3.02 1.70 1.74 1.19 1.27 
草地 22.86 21.34 12.49 11.52 10.37 9.82 
城乡 工矿 居民 用 地 1.95 2.40 1.14 1.40 0.80 1.00 
水 域 5.33 3.79 4.07 2.45 1.28 1.35 
未 利用 土地 94.33 95.10 80.61 81.64 13.72 13.46 
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(a) 土地 利用 变化 
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(c) 阿 姆 河流 域 及 分 区 年 蒸 散 量 
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图 7 成 海流 域 土地 利用 变化 百分比 和 蔡 散 量变 化 


Fig.7 Change percentage of land use and change of evapotranspiration in the Aral Sea Basin 
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区 域 水 储量 减少 ;而 以 漫灌 为 主 的 灌溉 方式 也 可 能 
市 来 灌区 地 下 水 位 的 抬升 ,使 区 域 水 储量 增加 。 但 
总 体 而 言 ,耕地 面积 扩张 带 来 的 水 资源 耗 散 增 大 ， 
陆地 水 储量 以 减 小 为 主 ; 在 局 部 区 域 ,如 阿 姆 河流 
域 中 游 地 区 ,地 下 水 位 抬升 和 水 库 蓄 水 增加 可 能 是 
导致 区 域 陆 地 水 储量 增加 的 主要 原因 。 流 域 尺 度 
上 ,流域 的 大 部 分 地 区 蒸 散 量 呈 增加 的 趋势 ,其 中 
整个 锡 尔 河流 域 和 阿 姆 河流 域 的 中 游 获 散 量 大 致 
上 时 增 加 趋势 ; 阿 姆 河 流域 上 游 山 区 及 下 游 的 西部 
蒸 散 量 呈 减少 趋势 (图 ga) 。 蒜 散 量 与 陆地 水 储量 
变化 的 相关 性 强 于 降水 量 和 气温 ,高度 相关 的 区 域 
在 阿 姆 河流 域 中 下 游 交 界 处 的 灌区 ;人 负 相 关 地 区 在 


(a) RE LH 


威海 流域 下 游 西 南部 和 锡 尔 河流 域 中 下 游 交 界 处 
(图 8b) , 茸 散 量 的 增加 是 土地 利用 变化 和 区 域 升温 
共同 影响 的 结果 。2002 一 2016 年 各 要 素 变 化 率 反 
映 了 威海 流域 陆地 水 储量 呈 减 少 态 势 不 是 单一 因 
素 造 成 的 。 降 水 量 减少 ,气温 升 高 导致 蒸 散 增加 ， 
灌溉 农业 的 发 展 、 耕 地 面积 增加 导致 水 资源 耗 散 ， 
综合 因素 导致 威海 流域 整体 的 陆地 水 储量 呈 减 少 
趋势 ( 表 2)。 


4 结论 与 讨论 


本 文 利 用 GRACE 数据 提取 咸 海流 域 陆地 水 储 
量变 化 时 间 序 列 ,并 结合 区 域 气候 要 素 和 土地 利用 
数据 ,探讨 了 气候 变化 与 人 类 活动 对 流域 水 储量 变 


(b) 蒸 散 量 与 陆地 水 储量 关系 


mu 0.50~0.74 


图 8 成 海流 域 磁 散 量 变化 及 其 与 陆地 水 储量 的 关系 
Fig. 8 Chang of evapotranspiration in the Aral Sea Basin and its relationship with TWS 
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#2 2000 一 2016 年 各 要 素 的 变化 率 


Tab. 2 Change rates of various elements from 2000 to 2016 


流域 EK Et /mm- a" 气温 /Ca 陆地 水 储量 /mm.a'  2XBLE/I0 m-a" ”耕地 /10' km? 
成 海流 域 -1.14 0.110 -3.20 21.63 1.65 
阿 姆 河 流域 上 游 -2.14 0.002 -7.28 3.80 0.04 
中 游 -1.07 0.010 1.79 8.48 0.76 
下 游 0.27 0.005 -7.36 -0.26 0.39 
锡 尔 河流 域 Lv -3.58 -0.003 -3.81 4.33 0.15 
中 游 -0.98 0.001 -6.48 2.11 -0.01 
下 游 -1.02 0.001 -5.44 2.84 0.31 
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化 的 影响 。 得 出 如 下 主要 结论 : 

(1) 2002 一 2016 年 威海 流域 陆地 水 储量 总 体 呈 
亏损 状态 ,其 中 流域 西北 部 威海 地 区 陆地 水 储量 亏 
员 最 为 严重 , 主要 原因 是 办 联 解体 后 乌 效 别克 斯 
坦 继续 强化 农业 经 济 , 阿 姆 河中 下 游 农业 灌溉 面 

让 继续 扩张 ,加 之 持续 升温 的 影响 ,区 域 蒸 散 量 耗 
水 增加 “ 。 

(2) 陆地 水 储量 变化 与 降水 量 的 相关 系数 高 于 
地 表 温 度 ,气温 与 陆地 水 储量 的 相关 性 在 流域 大 部 
分 地 区 并 不 显著 ,说 明 地 表 温 度 不 是 陆地 水 储量 变 
化 的 主导 因素 。 

(3) 2000—2015 年 威海 流域 内 农作物 种 植 面 积 
增加 和 灌溉 导致 水 资源 耗 散 增加 ,加剧 了 成 海流 域 
水 量 输出 。 戌 海流 域 年 燕 散 量 以 21.63x10' m'a 3E 
率 增 加 , 且 与 陆地 水 储量 的 变化 相关 性 较 强 (最 大 
相关 系数 0.74) ,是 影响 陆地 水 储量 变化 的 主要 原 
之 一 ;但 由 于 蒸 散 量变 化 的 驱动 要 素 不 仅 包 括 气 
温 、 风 速 .相对 湿度 等 气象 因素 ,也 与 农作物 .天然 
植被 及 水 体 的 面积 和 时 空 分 布 密切 相关 ,影响 机 制 
复杂 ,在 目前 的 数据 基础 上 很 难 定量 区 分 蒸 散 量变 
化 是 否 是 陆地 水 储量 的 主要 影响 因素 。 

(4) 南 咸 海 地 区 由 于 水 面 持续 萎缩 ,湖泊 水 储 
量 不 断 减 小 ,地 下 水 在 补充 水 体 蒸发 的 同时 也 被 抽 
取 用 于 农业 灌溉 ,导致 阿 姆 河 下 游 地 区 陆地 水 储量 
较 大 幅度 减 小 ; 而 北 威海 一 侧 ,湖泊 面积 和 水 量 均 
小 幅 增 加 ,但 水 体 蒸 发 也 随 湖面 扩张 和 升温 而 增 
加 ,加 之 农田 天 然 植被 和 湿地 蒸 散 量 耗 水 增加 , 
尔 河流 域 下 游 地 区 陆地 水 储量 呈 减 小 趋势 ,但 减 小 
幅度 低 于 阿 姆 河 流域 下 游 。 虽 然 本 研究 采用 的 是 
较 新 版 本 的 MASCON 解数 据 ,但 其 对 咸 海流 域 陆地 
水 储量 的 反 演 依然 存在 不 确定 性 ;同时 ,单一 模型 
模拟 的 燕 散 量 数据 不 确定 性 也 较 大 ,在 后 续 人 研究 


中 ,将 综合 RL06 和 RL05 多 版 本 陆地 水 储量 数据 和 
GLDAS 中 多 模型 模拟 的 蒸 散 发 数据 ,开展 数据 比较 
与 分 析 。 
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Abstract: The Aral Sea, which was once the largest inland lake in Central Asia, is located far from oceans, and 
thus, a vulnerable hydrological system. The water balance in the lake has changed in recent decades under the 
combined impacts of climate change and anthropological activities, and problems related to limited water 
resources and ecological environment have attracted wide attention. Analysis of the variation trend and 
influencing factors of terrestrial water storage (TWS), one of the most important parameters of land water 
variation, can provide a scientific basis for the sustainable utilization of water resources in the Aral Sea Basin. In 
this paper, the spatiotemporal variations of terrestrial water storage change (TWSC) in the Aral Sea Basin are 
evaluated using the JPL- RLO6M data collected by the GRACE gravity satellite from 2002 to 2016. CRU TS4.03 
meteorological data, GLDAS-Noah surface evapotranspiration data, and high-precision land use data are combined 
to assess time-series changes in TWSC, precipitation, temperature, and water evaporation via the linear trend 
method. Using a resolution of 1? x 1? pixel, we calculated the correlation coefficients of various parameters, 
including TWSC, precipitation, and surface temperature, to explore the effects of climate change and 
anthropological activities on TWS in the Aral Sea Basin. Several interesting results were obtained. (1) The TWS 
of the Aral Sea Basin declined at a rate of -3.20 mm + a`' in 2002—2016. A water surplus in spring and summer 
and a water deficit in autumn and winter were observed. Moreover, the TWSC of the Aral Sea Basin was 
characterized by a surplus in the central and eastern regions and deficits in peripheral areas. (2) Precipitation in 
the Aral Sea Basin in 2002—2016 decreased at a rate of — 1.14 mm-a , and the surface temperature increased at a 
rate of 0.11 ?C-a '. Precipitation and TWSC were positively correlated in most areas of the Aral Sea Basin but 
negatively correlated in the middle and lower reaches of the Amu Darya Basin and the lower reaches of the Syr 
Darya Basin. The surface temperature was negatively correlated with TWSC in most areas, but a positive 
correlation between these parameters was observed in some areas in the lower reaches of the Amu Darya and Syr 
Darya River Basins. TWSC was more strongly correlated with precipitation than with surface temperature. (3) In 
2000—2015, the cultivated land area in the Aral Sea Basin increased slightly by 1.65x 10* km?, but the grassland 
and water areas decreased. Increases in water for crop management and irrigation intensified water expenditures 
from the Aral Sea Basin. Under the joint influence of climate change and anthropological activities, 
evapotranspiration in the Aral Sea Basin increased at a rate of 21.63x10* m .a '. The correlation coefficient between 
evapotranspiration and TWSC was as high as 0.74, which means this parameter may be a main factor affecting 
TWSC. 
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